Charged bottomonium-like states Z fo (10610) and Z b ( 10650) and the T(55) T(2S)n + n- decay 



O 

(N 



43 
i , 

Ph. 
43 



> 

m 

OS 

in 

in 
o 



X 



Dian-Yong Chen u EI Xiang Liu 12t | and Shi-Lin Zhu 4 fj 

1 Research Center for Hadron and CSR Physics, Lanzhou University and Institute of Modern Physics of CAS, Lanzhou 730000, China 
2 School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China 
3 'Nuclear Theory Group, Institute of Modern Physics of CAS, Lanzhou 730000, China 
4 Department of Physics and State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing 100871, China 

(Dated: January 25, 2013) 

Inspired by the newly observed two charged bottomonium-like states, we consider the possible contribution 
from the intermediate Z/,(10610) and Z>,(10650) states to the T(5S) — > Y(2S)n + n~ decay process, which nat- 
urally explains Belle's previous observation of the anomalous Y(2S )n + n~ production near the peak of T(55) 
at = 10.87 GeV [K.F. Chen et al. (Belle Collaboration), Phys. Rev. Lett. 100, 112001 (2008)]. The re- 
sulting dT(T(55 ) -» y{2S )7i*n-)ldm„ + „- andc/T(Y(5S) -» Y(2S>r + 7r)/d cos distributions agree with Belle's 
measurement after inclusion of these Z& states. This formalism also reproduces the Belle observation of the 
double-peak structure and its reflection in the Y(2S)7r + invariant mass spectrum of the Y(5S) — > Y(2S)?r + 7r~ 
decay. 



PACS numbers: 14.40.Pq, 13.25.Gv, 13.66.Bc 



Very recently, the Belle Collaboration announced the 
first observation of two charged bottomonium-like states 
Z/,(10610) and Z/,(10650) in the hidden-bottom decay chan- 
nels Y(«5)tt ± (n = 1, 2, 3) and h^mP)^ (m = 1, 2) of Y(5S) 
[Q]]. The measured parameters of Z/,(10610) and Z/,(10650) 



M Zij(1 o6io)/r Zi ,(i06io) = 10608.4+ 2.0/15.6 ± 2.5 MeV, 
M Zi(1 0650)/r Zi ,(i0650) = 10653.2+ 1.5/14.4+ 3.2 MeV. 

The analysis of the angular distribution indicates that the 
quantum numbers of both Z/,(10610) and Zft(10650) are 
lG(J p ) = 1+(1+). Both Z fe (10610) and Z&(10650) are charged 
hidden-bottom states. Moreover they are very close to the 
thresholds of BB* and B*B* (H, respectively. Thus, Z fe (10610) 
and Zfc(10650) are ideal candidates of the BB* and B*B* S- 
wave molecular states, which were studied extensively in 
Refs. SI. 

On the other hand, a new puzzle arises in the theoreti- 
cal study JH |a] of the dipion invariant mass distribution and 
the cos 6 distribution of the anomalous T{2S produc- 
tion near the peak of Y(5S) OD. While all the other cal- 
culations are well in accord with the Belle data, the pre- 
dicted differential width dT{f{5S) -> Y(2S)7r + 7r)/dcos6> 
disagrees with the Belle measurement JH. In this work, we 
will illustrate that the inclusion of these two Zb states in the 
Y(55 ) — > Y(2S )n + n~ decays explains the puzzling line shape 
of dF(r(5S ) -> T(2S )n + n-)/d cos 6 very naturally. 

In general, there exist three mechanisms for the Y(5S) 
hidden-bottom decays with the dipion emission 

Y(5S)^ r{2S)(p l )n + {p 1 )n-{ Pi ). 
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The first one is the Y(2S);r + ;r~ direct production by Y(5S) 
decay (see Fig. Q] (a)). The so-called direct production 
of Y(5S) — » Y(2S )7r + 7r~ denotes that there does not ex- 
ist the contribution from the intermediate mesons (such as 
o-(600),/ (980), hadronic loop constructed by B {t) or B ( * } 
mesons, Zb) to Y(5S) — > T{2S)n + ji~ . Thus, the direct produc- 
tion of Y(5S ) — > T{2S provides the background contri- 
bution. 

The QCD Multipole Expansion (QME) method JH is gen- 
erally applied to deal with the dipion transition between heavy 
quarkonia. So far, there exist many theoretical efforts study 
the dipion transitions between the bottomonia ITlUl8ll (see 
Refs. 1IT9I - I2TI1 for a detailed review). In this work, we do 
not intend to calculate the contribution from the direct tran- 
sition under the framework of the QME method, but alter- 
natively follow the effective Lagrangian approach to describe 
Y(5S) — > T{2S)n + ji~ transitions. The transition amplitude of 
the direct production of Y(5S ) — > Y(2S )7i + n~ can be written 
as 



Wt[Y(5S) -> T(2S)7!- + 7T] 

Direct 
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which was suggested by Novikov and Shifman in the study 
of the if/' — > J/i[m + n~ decay I22I1 . where the subscripts S- 
wave and D-wave denote the S-wave and D-wave contribu- 
tions respectively. AM is the mass difference between T(5S ) 
and T(2S). q 2 = (p2 + Pi) 2 = m 2 +ji _ denotes the invariant 
mass of 7r + 7r~, while 8 is the angle between T(5S) and tC in 
the 7T + 7r~ rest frame. The pion decay constant and mass are 
taken as f„ - 130 MeV and m n = 140 MeV, respectively. In 
Eq. ([TJ, K and T are free parameters to be determined when 
fitting the experimental data. 

Different from the other low-lying bottomonia with J PC = 
1 , T(5S ) is above the B ( * , B ( * > thresholds and predominantly 
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FIG. 1: The diagrams in the T(55) hidden-bottom decay. Here, Fig. 
[2(a) represents the T(55) direct decay into Y(2S )n + n~ , while Fig.[T] 
(b) denotes the intermediate hadronic loop contribution to T(55 ) — » 
T(25 )n + n~ . (c) and (d) describe the intermediate Z ; * contribution to 
Y(55) -> r(2S)n + n-, where Z± = (Z^lOeiO)*, 26(10650)*}. 



decays into flWfiW p a j ri which may render the coupled chan- 
nel effect quite important fjl-Holl. When exploring the Y(55') 
hidden-bottom decay, the coupled channel effect has to be 
taken into account. In other words, there also exists the sec- 
ond mechanism contributing to the T{5S) — > T{2S)n + n~ 
transitions as shown in Fig. Q] (b), where the intermediate 
and hadronic loop is the bridge to connect the ini- 
tial state T(5S) and final state T(2S)n + n~. Furthermore, 
Y(5S) — » Y(nS)ji + n~ can be approximately expressed as a 
sequential decay process. Y(5S) first transits into Y(2S) and 
the scalar meson <x(600). Then cr(600) couples with the dip- 
ion. Choosing cr(600) as the intermediate state contribution to 
the Y(5S) — > Y(2S)7r + 7r~ process is not only consistent with 
the Belle data 13,01 but also allowed by the phase space of the 
decay channel. 

If comparing the dipion invariant mass spectrum of 
T(5S) -> Y(2S)7r + 7T in Refs. HH, the data in Ref. H at 
the higher end of m„+„- are qualitatively different from those 
in Ref. 0], where the total events in Ref. J3l are at least one 
order of magnitude less than those in Ref. W. Such a large 
accumulation of events at m^ K - > 700 MeV Ol might be due 
to the contribution from the tail of the intermediate /o(980). 
We did not include the /o(980) contribution when we analyzed 
the data in |7||. Considering the situation of the new data of 
the dipion invariant mass spectrum of Y(55 ) — > Y(25 )n + n~ 
[Q]], we also include /o(980) contribution to the analysis of 
Y(55) -> Y(2S)7r + 7T in the following. 

The effective Lagrangians relevant to the Fig. [2(b) include 

-Crss = igresV»(d f '£S f - BPS*), (2) 
£rs-s = -ito•8^"" i 3„Tv(W;S t +S^ , ), (3) 

+(<9 / ,y v S* v - Y v 5 jU S* v )S vt 

+S v (T v^ s .t _ ^ TyS *vt )}j ( 4 ) 



and 



(5) 

where S = (B°, B , B~) and (&) T = (B°, B + , B+). There are 4 
diagrams. Thus, the concrete expressions of decay amplitudes 



are written as 



■^-^[igr(5S)BB4(5S)^P^ - iPiJH 

x[igr(nS)BB^. (nS) (-ipi P - iq P )][gBBs] 
1 1 1 



2 2 _ 2 2 
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J (2n) A ^~ 8r(SS >BB ' e wp(~iPc) e r(5S ) (iPi)] 



x[-gr{nS)BBeg T 94,(ip 6 3 )e^ (nS) (iq e )][-gB'B's] 
1 ~g Pp + v\f % lm\ -g* p + q*q p /m 2 B , 
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^-^[-igT(5S)B-B'4 (5s) ((ip2^ ~ ipi M )g Vf 

+(-ipo P - ip2p)g flv + (ipiv + ipov)gp P )] 

X[-igr(nS)B'B'^(nS)^~'P l <^ ~ ^)gaB 

+(ip3B + ip\B)ga^, + (iq a - ipia)gB<t,)][-gB'B's] 

,-g pa + f x PV<- s VT + pIpIK 



X- 
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P2 

-g BT + q B q T lml, 
q z - m B , 



(9) 



The amplitude M C AB indicates that the initial T(5S) dissolves 
into intermediate AB, which transit into the final Y(25 ) and 
scalar meson by exchanging meson C. In the above ex- 
pressions, the form factor is introduced by T{q 2 ) = (A 2 - 
m 2 E )/(q 2 - m 2 E ). And ni£ is the mass of the exchanged B (t) 
meson in the B i -* ) B { * ) — > Y(2S)S transitions shown in Fig. Q] 
(b) and A = m E + aA QC D with Aq C d = 220 MeV. As indi- 
cated in Ref. Hf], we can parameterize the decay amplitude of 
T(5S) — > Y(2S)7i + n~ corresponding to Fig. Q](b) as 



M[Y(5S) ^ Y(25)cr(600) ^ Y(25)^ + ^ _ ] 
e* } 



Ciy>/T (■ \ F (J 



(P2 + Pi) 2 ~ ml + imo-Ttj- ' 



(10) 



if only considering the S-wave contribution. Here, we intro- 
duce F a as the fitting parameter. 

Similar to Eq. (llOl l. the parameterized decay amplitude of 
Y(5S) T{2S)n + n- with / (980) as the intermediate state 
can be expressed as 



M[T(5S) -» Y(2S)/ (980) -» Y(25)^ + tt- 

e T (5S) • 4(25)^/0 

(P2 + /5 3 ) 2 - m 2 + im fo r fo ' 



(ID 
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which corresponds to Fig. [T](b) with the replacement <x — > 
/o(980). 

Regarding the contribution of these two newly observed 
Zb states to the Y(5S) — > Y(2S )7r + 7r~ process, we introduce 
the third mechanism depicted in Fig. Q] (c) and (d), where 
Z^s are the intermediate states and interact with T(5S)7r T and 
T(25')^ ± . The general expressions of the amplitudes of Fig. [TJ 
(c) and (d) are 



z b T(5S) e T(2S) 



-g^v + (K+^)(Pl+PP/< 

(/?i + Pi) 2 - m 2 z + !'m Z6 r Z( , 



M[T(5S) 



Z^tt 



= F zr e£ 



v7 c T(5S) c T(2S)" 



T(25)7T^ + ] Z - 
Q?l + P3) 2 - m| + !'m Z6 r Z(; 



(12) 



,(13) 



and 



respectively, where we define F z + = g T(5S)Z+ „,? z+T(2S> , + 
F Z7, = S y{$S)Z - b „g z - yaS) „- ■ Since Fi §- 03 ( c ) and ( d ) ^ related 
to each other by charge-conjugation, thus F z - = F z + = Fz b ■ 

Thus, the total decay amplitude of the Y(5S ) — > Y(2S )7r + 7r~ 
decay is 

Direct 

+e*M[T(55) -> T(25)cr(600) Y(2S)7r + 7T] 

+e^/oM[T(55) -> T(25)/ (980) -> Y(2S KtT] 
+ J] e*^ (m[T(55) -» Z+tT -» Y(25)7r + 7r-] Z): 



+M[T(55) -> Z,;tt + T(25)7r + 7r-] z - 



(14) 



where we have introduced the phase angles (pa-, <pf a , <fi Zb(106l0) 

aIld ^,,10650,- 

As a three body decay, the differential decay width for 
Y(5S) T(25)tt + ^- read as, 



i 



3 (2tt) 3 32m3, (5s) 



|At tot ail 2 ^ (2S)?r dm^ (15) 



with m 2 (25 . )7r+ = (p\+p2) 2 and;7z^ +jr _ = {pi+pi) 2 - The relevant 
resonance parameters are listed in Table. [TJ 

TABLE I: The resonance parameters adopted in our calculation U, 



State 


Mass (GeV) 


State 


Mass (GeV) Width (GcV) 


Y(55) 


10.870 


o-(600) 


0.478 0.324 






/o(980) 


0.980 0.100 


Y(25) 


10.023 


Z,,(10610) 


10.608 0.0156 






Z b ( 10650) 


10.653 0.0144 



If considering only the contributions from Fig. [2(a) and 
(b) in our present scenario, we have four free parameters as 
listed in Table [TTJ where the <t(600) contribution is included 
to fit the Belle data 0]. With the help of the MINUIT pack- 
age, we perform the global fit to the experimental data of the 



dipion invariant mass spectrum distribution and the cos 9 dis- 
tribution of the Y(2S production near the peak of Y(5>S ) 
0]. The best fit to the dipion invariant mass spectrum dis- 
tribution is shown in the left-panel in Fig. [2] Unfortunately 
the corresponding cos 6 distribution of the Y(2S )n + n~ pro- 
duction strongly deviates from the Belle data as shown in the 
right-panel of Fig. [2] The values of the obtained fitting pa- 
rameters are presented in Table [EI] Such discrepancy between 
theoretical and experimental results stimulates a New Puzzle 
first indicated in Ref. JH. At present, solving these new 
puzzle becomes an important and intriguing research topic, 
which will be helpful to underlying mechanism behind the 
Y(5S) Y(2S )n + n- decay. 
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FIG. 2: (Color online). The dipion invariant mass (m nn ) distribution 
(left-panel) and the cos# distribution (right-panel) of the T(55) — » 
T(25 )n + n~ decay. The dots with error bars are the results measured 
by Belle 0], while the green histograms is the best fit from our model 
without including the intermediate Z;,(10610) ± and Zf,(10650) ± con- 
tribution to Y(5S) — > Y(2S)7r + ;r~. When fitting the experimental 
data 1 7], we only include cr(600) contribution. 



TABLE II: The values of the fitting parameters for the best fit to the 
Belle data of Y(25 production near the peak of T(55 ) 0] with- 
out considering the contributions from Z/j(10610) and Z/,(10650). For 
the obtained central values of the parameters, the corresponding par- 
tial decay width of Y(5S ) Y(2S )7r + 7T is 0.836 MeV. 



Parameter 



Value 



Parameter 



Value 



T 0.943 ±0.071 

25.603 ±2.175 GeV 2 



0.739 ± 0.034 
2.623 ± 0.132 Rad 



In contrast, we consider the contribution from Z/,(10610) 
and Z/,(10650) in the following and discuss the dependence 
of dT/dm„+„- and dF/dcosd of T(5S) Y(2S)7r + 7r on 
m„+„- and cos respectively. Under this scheme, we refit the 
Belle data 0}] with Eq. tl4) . There are 10 fitting parame- 
ters as listed in Table [HI] In Fig. [3j we present a compar- 
ison between the Belle data (dots with error bars) and our 
best fit (histograms) to the Belle data [1], which indicates 
that the line shapes of the invariant mass spectra of 7r + 7r~ and 
Y(2S)7r + for Y(5S) -> Y(2S)7r + 7r describe the Belle data 
10] well. The double-peak structure around 10.6 GeV and 
its reflection around 10.25 GeV are reproduced by our model 
well. With the central values of these parameters in Table UTTJ 
we obtain the partial decay width of Y(55) — > Y(2S)7t + 7t~ 
F = 0.915 MeV, which is consistent with the Belle measure- 
ment T = 0.85 ± 0.07(stat) ± 0.16(syst) MeV 0]. Thus, the 
contribution from these charged Zt, resonances provides a pos- 
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sible solution to the puzzle why the Y(5S ) — * T(2S )n + n de- 
cay width is abnormally large Q. 



TABLE III: The values of the fitting parameters for the Y(5S) — » 
Y(2S)jr + ;r~ decay after including the contributions from Z;,(10610) 
and Z b ( 10650). 



Parameter 


Value 


Parameter 


Value 


r 


1.404 ±0.068 


K 


0.301 ±0.013 


Fa 


20.037 ± 0.423 GeV 2 


(pa 


0.907 ± 0.132 rad 


Ff 


17.076 ± 3.563 GeV 2 


<Pfo 


-0.753 ± 0.140 rad 


^Z t (10610) 


3.412 ± 0.385 GeV 2 


*rz fc (10610) 


-3.135 ±0.030 rad 


fzj(10650) 


2.994 ± 0.261 GeV 2 




-2.836 ±0.165 rad 
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FIG. 3: (Color online). The invariant mass spectra of n + n~ and 
Y(25> + for T(55) -> Y(2S)7r + 7T. Here, the histograms are the- 
oretical results obtained in our scenario including the intermediate 
/o(980) contribution, while dots with error bars are the Belle data in 
Ref. Q|]. We also plot the / (980) contribution separately (red solid 
lines). 

From Table Hill we notice that the uncertainty of Ff is one 
order-of-magnitude larger than that of F a , which means the 
fit is less sensitive to the /o(980) than to the <x(600). Using 
Eq. ( fT4l . we reanalyze the new Belle data in Ref. with the 
obtained fitting parameters in Table [IV] where we do not in- 
clude the /o(980) contribution. The comparison between our 
fitting result and the experimental data are given in Fig. |4] By 
the scenario in Eq. (fT4l . we reproduce the Belle data well, 
which confirms that the intermediate /o(980) contribution to 
T{5S) — > Y(2S)n + n~ is small. If comparing the obtained 
values of the fitting parameter in Tables UTU and |IV] we no- 
tice that the eight common parameters do not change much in 
the two schemes. With the parameters listed in Tables HIT1 and 
IIV1 we also present the cos 9 distribution with and without the 
intermediate /o contribution. The experimental measurement 
of the cos 9 distribution for Y(5S ) — > T(2S )n + n~ J7[| can be 
described well with the scenarios in this work. This fact in- 
dicates that the two Z/, structures play important role in the 
understanding of the Belle data, especially the cos 9 distribu- 
tion of Y(55) -> T(2S)n + 7T-. 

In summary, the Belle Collaboration announced an ex- 
citing observation of two charged bottomonium-like states 
Z/,(10610) and Z\,( 10650). These Z/, states are good candidates 
of exotic states, which calls for theoretical efforts in revealing 
their underlying structures. Carrying out the phenomenolog- 
ical study relevant to Z/,(10610) and Z/,(10650) is one of the 
important and valuable issues of heavy quarkonium physics, 
which is full of challenges and opportunities l24ll25ll . 



TABLE IV: The values of the fitting parameters for the Y(5S) — » 
T(2S)n + n~ decay after including the contributions from Z;,(10610) 
and Z b { 10650). These parameters are obtained by fitting the new 
experimental data without including the contributions from /o(980). 



Parameter 


Value 


Parameter 


Value 


T 


1.073 ±0.064 


K 


0.379 ± 0.034 


Fa 


23.833 ± 2.503 GeV 2 


<Pa 


1.127 ±0.128 rad 


Fz b (W6W) 


3.200 ± 0.345 GeV 2 


^Z fc (10610) 


-3.141 ±0.076 rad 


-Fz 6 (10650) 


2.686 ± 0.306 GeV 2 


^Z,,( 10650] 


-2.703 ± 0.225 rad 
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FIG. 4: (Color online). The distribution invariant mass spectra m n „ 
and m raS) „ for Y(55) — > Y(2S)7r + ;r~ without including the contribu- 
tions from / (980). Here, we use Eq. dl4t to redo the analysis. The 
histograms are the fitting results. The dots with errors correspond to 
the Belle data (J. 



The Zi(10610) and Z/,(10650) states are related to the 
anomalous phenomena of Y(2S)n + n~ production near Y(55 ) 
previously reported by Belle [7]. Comparing the fitting re- 
sults without and with the contributions from the newly ob- 
served states, we notice that the intermediate Z/,(10610) and 
Zi,( 10650) play a crucial role in the behavior of c/F(Y(55) — > 
Y(2S)7r + 7r)/flfcos0. The inclusion of the Z/,(10610) and 
Z fe (10650) contribution to Y(5S) Y(25)tt + ^ provides a 
unique mechanism of understand the puzzling cos 9 distribu- 
tion of Y(2S )7r + 7r~ production near T(5S) 0]. The double- 
peak structure and its reflection in the Y(2S invariant mass 
spectrum of Y(55) — > Y(2S)ji + tt~ 0]] are also reproduced by 
this mechanism. In this work, the values of the fitting param- 
eters in our scenario are obtained by fitting Belle data [OJ |7| ■ 
To some extent, the interpretation of the values of these pa- 
rameters is related to the understanding of background, the 
structures of two Zb states etc, which is an interesting research 




cos(0) cos(0) 



FIG. 5: (Color online). The cos0 distributions for Y(55) — » 
Y{2S)7i + 7i~ . The histograms in the left-hand side and right-hand side 
diagrams are the fitting results without and with the contribution from 
/o(980). The dots with errors correspond to the Belle data @). 
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topic. 

Besides finding the signals of Zj,(10610) and Z/,(10650) 
in T(25')^ ± decay channel, Belle's analysis of its remain- 
ing four hidden-bottom decay channels Y(nS )n ± (n = 1 , 3) 
and hbimPJn* (m = 1,2) also indicate the observation of 
Z/,(10610) and Zft(10650) The present formalism can 
be extended to study the dipion invariant mass distribution 
and the cos 6 distribution of Y(5S) — > Y(15\ ?>S)ji + ji~ and 
T(55) -> h b (\P,2P)n + n- decay. 

Additionally, Belle's measurement favors the BB* and B*B* 
molecular explanation of the Z/,(10610) and Zj,(10650) res- 
onances respectively. The possible S-wave BB* and B*B* 
molecular states were investigated extensively in Refs. JH3]. 
Very recently, the authors in Ref. | j26j] discussed the special 
decay behaviour of the J=l S-wave BB* and B*B* molecular 
states based on the heavy quark symmetry. Future dynamical 



study of the mass and decay pattern of the S-wave BB* and 
B*B* molecular states are very desirable. 
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